(* Content-type: application/vnd.wolfram.mathematica *)

(*** Wolfram Notebook File ***)
(* http://www.wolfram.com/nb *)

(* CreatedBy='Mathematica 10.0' *)

(*CacheID: 234*)
(* Internal cache information:
NotebookFileLineBreakTest
NotebookFileLineBreakTest
NotebookDataPosition[       158,          7]
NotebookDataLength[     38162,        922]
NotebookOptionsPosition[     36915,        873]
NotebookOutlinePosition[     37274,        889]
CellTagsIndexPosition[     37231,        886]
WindowFrame->Normal*)

(* Beginning of Notebook Content *)
Notebook[{

Cell[CellGroupData[{
Cell["Numeric Integration Test", "Title",
 CellChangeTimes->{{3.636142598084065*^9, 3.636142607964061*^9}}],

Cell[BoxData[{
 RowBox[{
  RowBox[{"SetDirectory", "[", 
   RowBox[{"NotebookDirectory", "[", "]"}], "]"}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"Porder", "=", "50"}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"workpre", "=", "30"}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"method", "=", 
   RowBox[{"{", 
    RowBox[{"\"\<GlobalAdaptive\>\"", ",", 
     RowBox[{"Method", "\[Rule]", "\"\<GaussKronrodRule\>\""}], ",", 
     RowBox[{"\"\<SymbolicProcessing\>\"", "\[Rule]", "0"}]}], "}"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"UseExactIntegration", "=", "False"}], ";"}]}], "Input",
 CellChangeTimes->{{3.636142615306436*^9, 3.636142647856892*^9}, 
   3.636142767163019*^9, 3.636142962846133*^9, 3.636142992896546*^9, {
   3.636143508128821*^9, 3.636143521421315*^9}, {3.63614454650169*^9, 
   3.636144546794828*^9}, {3.636144901934841*^9, 3.63614490422192*^9}, {
   3.636145660979344*^9, 3.636145668485813*^9}, {3.6361457095649776`*^9, 
   3.636145712964238*^9}, {3.636145763563759*^9, 3.636145779745656*^9}, {
   3.6361458145332413`*^9, 3.636145848735752*^9}, {3.636145937889634*^9, 
   3.636145946495983*^9}, {3.636146167417938*^9, 3.636146174434998*^9}, {
   3.636146270587102*^9, 3.6361462708765287`*^9}, {3.636146419886677*^9, 
   3.636146436165842*^9}, {3.636146545151203*^9, 3.63614654522635*^9}, {
   3.636148158037044*^9, 3.63614816979571*^9}, {3.636148234490814*^9, 
   3.6361482346024923`*^9}, {3.636154994984722*^9, 3.6361550083751593`*^9}, {
   3.6361551317279*^9, 3.636155131973566*^9}, {3.636155535336226*^9, 
   3.636155592160096*^9}, {3.636155631582512*^9, 3.636155643114778*^9}, {
   3.636155695961226*^9, 3.636155707544558*^9}, 3.63629367693571*^9, {
   3.6362940048974657`*^9, 3.636294005019759*^9}, 3.636294098895577*^9, {
   3.636294145296986*^9, 3.636294145388741*^9}, {3.636294517952807*^9, 
   3.636294518047888*^9}, {3.636294729175123*^9, 3.636294729250318*^9}, {
   3.6362950607639933`*^9, 3.636295060871091*^9}, {3.636295125079165*^9, 
   3.636295154326345*^9}, {3.636295404107905*^9, 3.636295404826199*^9}, {
   3.636295731405182*^9, 3.636295734236137*^9}, {3.636296015176161*^9, 
   3.636296019579637*^9}, {3.6362965329274178`*^9, 3.636296534703166*^9}, 
   3.6363120028577538`*^9, {3.636312496694271*^9, 3.63631250416012*^9}, {
   3.636312536788467*^9, 3.636312536898303*^9}, 3.636312627820691*^9, {
   3.636314284068304*^9, 3.636314308730073*^9}, {3.636316391043355*^9, 
   3.636316394081561*^9}, {3.6363197795256166`*^9, 3.636319785228339*^9}, 
   3.636391463648078*^9, {3.636392533390601*^9, 3.636392533497951*^9}, 
   3.636393155024577*^9, {3.636393208511507*^9, 3.636393208590991*^9}, {
   3.636393550228066*^9, 3.6363935503377247`*^9}, {3.636393819128875*^9, 
   3.6363938192466087`*^9}, {3.636398000516883*^9, 3.636398002076282*^9}, 
   3.636398032689333*^9, {3.636398373611994*^9, 3.636398373747633*^9}, {
   3.636398498779296*^9, 3.636398498872786*^9}, {3.636398812587636*^9, 
   3.636398812649754*^9}}],

Cell[CellGroupData[{

Cell["Functions", "Section",
 CellChangeTimes->{{3.636295532894163*^9, 3.6362955354443073`*^9}}],

Cell[BoxData[{
 RowBox[{
  RowBox[{"onedimensional", "=", 
   RowBox[{
    RowBox[{"(", 
     RowBox[{"i", "+", "1"}], ")"}], " ", 
    SuperscriptBox[
     RowBox[{"(", 
      RowBox[{"\[Xi]", "+", 
       FractionBox["1", "10"]}], ")"}], "i"]}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"twodimensional", "=", 
   RowBox[{"If", "[", 
    RowBox[{
     RowBox[{
      RowBox[{"i", "+", "j"}], "\[LessEqual]", " ", "P"}], ",", 
     RowBox[{
      RowBox[{"(", 
       RowBox[{"i", "+", "1"}], ")"}], " ", 
      SuperscriptBox[
       RowBox[{"(", 
        RowBox[{"\[Xi]", "+", 
         FractionBox["1", "10"]}], ")"}], "i"], " ", 
      RowBox[{"(", 
       RowBox[{"j", "+", "1"}], ")"}], " ", 
      SuperscriptBox[
       RowBox[{"(", 
        RowBox[{"\[Eta]", "+", 
         FractionBox["1", "10"]}], ")"}], "j"]}], ",", "0"}], "]"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"threedimensional", "=", 
   RowBox[{"If", "[", 
    RowBox[{
     RowBox[{
      RowBox[{"i", "+", "j", "+", "k"}], "\[LessEqual]", " ", "P"}], ",", 
     RowBox[{
      RowBox[{"(", 
       RowBox[{"i", "+", "1"}], ")"}], " ", 
      SuperscriptBox[
       RowBox[{"(", 
        RowBox[{"\[Xi]", "+", 
         FractionBox["1", "10"]}], ")"}], "i"], " ", 
      RowBox[{"(", 
       RowBox[{"j", "+", "1"}], ")"}], " ", 
      SuperscriptBox[
       RowBox[{"(", 
        RowBox[{"\[Eta]", "+", 
         FractionBox["1", "10"]}], ")"}], "j"], " ", 
      RowBox[{"(", 
       RowBox[{"k", "+", "1"}], ")"}], " ", 
      SuperscriptBox[
       RowBox[{"(", 
        RowBox[{"\[Zeta]", "+", 
         FractionBox["1", "10"]}], ")"}], "k"]}], ",", "0"}], "]"}]}], 
  ";"}]}], "Input",
 CellChangeTimes->{{3.636295537330488*^9, 3.6362955867999*^9}, {
   3.63631693891726*^9, 3.6363169392359457`*^9}, {3.6363169843802633`*^9, 
   3.636317058198694*^9}, {3.6363171966231937`*^9, 3.636317250277629*^9}, {
   3.636391706438459*^9, 3.6363917149138927`*^9}, {3.6363917831713943`*^9, 
   3.636391804187422*^9}, 3.636391848999147*^9, {3.636391959724174*^9, 
   3.636391971349429*^9}, {3.636392393041389*^9, 3.636392405793995*^9}}]
}, Open  ]],

Cell[CellGroupData[{

Cell["Line parametric space", "Section",
 CellChangeTimes->{{3.636142521436479*^9, 3.636142524139866*^9}, {
  3.6361439875675163`*^9, 3.63614399019761*^9}, {3.636145131261984*^9, 
  3.636145131469054*^9}}],

Cell[BoxData[{
 RowBox[{
  RowBox[{
   RowBox[{"If", "[", 
    RowBox[{"UseExactIntegration", ",", "\[IndentingNewLine]", 
     RowBox[{
      RowBox[{"OneDimensional", "=", 
       RowBox[{"Table", "[", 
        RowBox[{
         RowBox[{"Integrate", "[", 
          RowBox[{
           RowBox[{"Total", "[", 
            RowBox[{"Flatten", "[", 
             RowBox[{"Table", "[", 
              RowBox[{"onedimensional", ",", 
               RowBox[{"{", 
                RowBox[{"i", ",", "0", ",", "p"}], "}"}]}], "]"}], "]"}], 
            "]"}], ",", 
           RowBox[{"{", 
            RowBox[{"\[Xi]", ",", 
             RowBox[{"-", "1.0"}], ",", "1.0"}], "}"}]}], "]"}], ",", 
         RowBox[{"{", 
          RowBox[{"p", ",", "0", ",", "Porder"}], "}"}]}], "]"}]}], ";"}], 
     "\[IndentingNewLine]", ",", 
     RowBox[{"OneDimensional", "=", 
      RowBox[{"Table", "[", 
       RowBox[{
        RowBox[{"NIntegrate", "[", 
         RowBox[{
          RowBox[{"Total", "[", 
           RowBox[{"Flatten", "[", 
            RowBox[{"Table", "[", 
             RowBox[{"onedimensional", ",", 
              RowBox[{"{", 
               RowBox[{"i", ",", "0", ",", "p"}], "}"}]}], "]"}], "]"}], 
           "]"}], ",", 
          RowBox[{"{", 
           RowBox[{"\[Xi]", ",", 
            RowBox[{"-", "1.0"}], ",", "1.0"}], "}"}], ",", 
          RowBox[{"Method", "\[Rule]", "method"}], ",", 
          RowBox[{"WorkingPrecision", "\[Rule]", "workpre"}]}], "]"}], ",", 
        RowBox[{"{", 
         RowBox[{"p", ",", "0", ",", "Porder"}], "}"}]}], "]"}]}]}], "]"}], 
   ";"}], "\[IndentingNewLine]"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"file1DL", "=", 
   RowBox[{"OpenWrite", "[", "\"\<Line.txt\>\"", "]"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"Do", "[", " ", 
   RowBox[{
    RowBox[{"Write", "[", 
     RowBox[{"file1DL", ",", 
      RowBox[{"OneDimensional", "[", 
       RowBox[{"[", "n", "]"}], "]"}]}], "]"}], ",", 
    RowBox[{"{", 
     RowBox[{"n", ",", "1", ",", "Porder"}], "}"}]}], "]"}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"Close", "[", "file1DL", "]"}], ";"}]}], "Input",
 CellChangeTimes->{{3.6361371506620207`*^9, 3.636137158025509*^9}, {
   3.636137188085631*^9, 3.6361372251130342`*^9}, {3.636137286833458*^9, 
   3.636137301169579*^9}, {3.636140265113166*^9, 3.6361402950367193`*^9}, {
   3.636140424428979*^9, 3.6361404358015823`*^9}, {3.636140469075234*^9, 
   3.636140478182131*^9}, {3.636140624330152*^9, 3.636140637423913*^9}, {
   3.6361410917859097`*^9, 3.636141131250258*^9}, {3.636141198603467*^9, 
   3.636141215139453*^9}, {3.636141249690695*^9, 3.6361412498709393`*^9}, {
   3.636141285406246*^9, 3.636141315181307*^9}, {3.6361416183823566`*^9, 
   3.6361416674598494`*^9}, {3.6361419265393553`*^9, 3.636141952905411*^9}, {
   3.636142010699959*^9, 3.6361420789895077`*^9}, {3.6361421105408297`*^9, 
   3.6361421950621758`*^9}, {3.6361422258123093`*^9, 
   3.6361423128862133`*^9}, {3.636142549876772*^9, 3.636142562616646*^9}, {
   3.636142690119768*^9, 3.636142754534297*^9}, {3.636142824855163*^9, 
   3.636142863735999*^9}, {3.6361429058220463`*^9, 3.636142958339533*^9}, {
   3.6361430242206373`*^9, 3.6361430279549017`*^9}, {3.6361431206765327`*^9, 
   3.636143264683196*^9}, 3.636143306946995*^9, {3.636143414419903*^9, 
   3.636143541667583*^9}, {3.63614361218032*^9, 3.636143705409419*^9}, {
   3.6361444077374067`*^9, 3.636144412445622*^9}, {3.636144538029191*^9, 
   3.636144543010495*^9}, {3.636144837112173*^9, 3.636144841098144*^9}, {
   3.636145788132073*^9, 3.636145801405088*^9}, {3.636145850706107*^9, 
   3.636145911201254*^9}, {3.636146073658416*^9, 3.63614611558902*^9}, {
   3.636146158424539*^9, 3.6361461787671003`*^9}, {3.636146297405676*^9, 
   3.636146297975884*^9}, {3.6361481833847647`*^9, 3.636148186070285*^9}, 
   3.636148223184607*^9, 3.6361556070064907`*^9, {3.636293663538445*^9, 
   3.636293668274581*^9}, {3.636294083480585*^9, 3.6362940893114347`*^9}, 
   3.636294136896103*^9, {3.636294973807578*^9, 3.636295017943521*^9}, {
   3.636295111974889*^9, 3.636295120714447*^9}, {3.6362951603964453`*^9, 
   3.6362951899052887`*^9}, {3.636295597744309*^9, 3.636295601617689*^9}, 
   3.636312164357861*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell["Triangle parametric space", "Section",
 CellChangeTimes->{{3.636142521436479*^9, 3.636142524139866*^9}, {
  3.636143971451137*^9, 3.6361439829008636`*^9}, {3.636144500980027*^9, 
  3.636144503019088*^9}, {3.6361451335567303`*^9, 3.636145137541813*^9}}],

Cell[BoxData[{
 RowBox[{
  RowBox[{"If", "[", 
   RowBox[{"UseExactIntegration", ",", "\[IndentingNewLine]", 
    RowBox[{
     RowBox[{"TwoDimensionalTriangle", "=", 
      RowBox[{"Table", "[", 
       RowBox[{
        RowBox[{"Integrate", "[", 
         RowBox[{
          RowBox[{"Total", "[", 
           RowBox[{"Flatten", "[", 
            RowBox[{"Table", "[", 
             RowBox[{
              RowBox[{"twodimensional", "/.", 
               RowBox[{"P", "\[Rule]", "p"}]}], ",", 
              RowBox[{"{", 
               RowBox[{"i", ",", "0", ",", "p"}], "}"}], ",", 
              RowBox[{"{", 
               RowBox[{"j", ",", "0", ",", "p"}], "}"}]}], "]"}], "]"}], 
           "]"}], ",", 
          RowBox[{"{", 
           RowBox[{"\[Xi]", ",", "0", ",", "1"}], "}"}], ",", 
          RowBox[{"{", 
           RowBox[{"\[Eta]", ",", "0", ",", 
            RowBox[{"1", "-", "\[Xi]"}]}], "}"}]}], "]"}], ",", 
        RowBox[{"{", 
         RowBox[{"p", ",", "0", ",", "Porder"}], "}"}]}], "]"}]}], ";"}], 
    "\[IndentingNewLine]", ",", "\[IndentingNewLine]", 
    RowBox[{
     RowBox[{"TwoDimensionalTriangle", "=", 
      RowBox[{"Table", "[", 
       RowBox[{
        RowBox[{"NIntegrate", "[", 
         RowBox[{
          RowBox[{"Total", "[", 
           RowBox[{"Flatten", "[", 
            RowBox[{"Table", "[", 
             RowBox[{
              RowBox[{"twodimensional", "/.", 
               RowBox[{"P", "\[Rule]", "p"}]}], ",", 
              RowBox[{"{", 
               RowBox[{"i", ",", "0", ",", "p"}], "}"}], ",", 
              RowBox[{"{", 
               RowBox[{"j", ",", "0", ",", "p"}], "}"}]}], "]"}], "]"}], 
           "]"}], ",", 
          RowBox[{"{", 
           RowBox[{"\[Xi]", ",", "0", ",", "1"}], "}"}], ",", 
          RowBox[{"{", 
           RowBox[{"\[Eta]", ",", "0", ",", 
            RowBox[{"1", "-", "\[Xi]"}]}], "}"}], ",", 
          RowBox[{"Method", "\[Rule]", "method"}], ",", 
          RowBox[{"WorkingPrecision", "\[Rule]", "workpre"}]}], "]"}], ",", 
        RowBox[{"{", 
         RowBox[{"p", ",", "0", ",", "Porder"}], "}"}]}], "]"}]}], ";"}]}], 
   "\[IndentingNewLine]", "]"}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"file2DT", "=", 
   RowBox[{"OpenWrite", "[", "\"\<Triangle.txt\>\"", "]"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"Do", "[", " ", 
   RowBox[{
    RowBox[{"Write", "[", 
     RowBox[{"file2DT", ",", 
      RowBox[{"TwoDimensionalTriangle", "[", 
       RowBox[{"[", "n", "]"}], "]"}]}], "]"}], ",", 
    RowBox[{"{", 
     RowBox[{"n", ",", "1", ",", "Porder"}], "}"}]}], "]"}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"Close", "[", "file2DT", "]"}], ";"}]}], "Input",
 CellChangeTimes->{{3.636143710281425*^9, 3.636143796641883*^9}, {
   3.6361438407822227`*^9, 3.636143866794297*^9}, {3.636143903500929*^9, 
   3.636143966659725*^9}, {3.63614407559863*^9, 3.6361440809457493`*^9}, 
   3.6361441283929443`*^9, {3.636144297998473*^9, 3.636144298326858*^9}, {
   3.636144377228163*^9, 3.636144447030261*^9}, {3.636144508874742*^9, 
   3.636144510005105*^9}, {3.636144566525465*^9, 3.636144612551203*^9}, {
   3.636144783860141*^9, 3.6361448083185167`*^9}, {3.636145903276352*^9, 
   3.63614590358461*^9}, {3.636146060562789*^9, 3.636146066882897*^9}, {
   3.636146099770914*^9, 3.636146116909381*^9}, 3.636146187887333*^9, 
   3.6361463025902452`*^9, {3.636148211364916*^9, 3.636148223155368*^9}, 
   3.6361556102685127`*^9, {3.636294180709875*^9, 3.6362941887225533`*^9}, {
   3.636294443843861*^9, 3.636294452828842*^9}, {3.636294745291925*^9, 
   3.636294750548738*^9}, {3.636294943877496*^9, 3.636294952537224*^9}, {
   3.636295034175301*^9, 3.6362950344145393`*^9}, {3.6362952154354*^9, 
   3.636295239615829*^9}, 3.6362953005090113`*^9, 3.636295392985839*^9, {
   3.636295505999547*^9, 3.63629550789007*^9}, {3.636295621479218*^9, 
   3.636295624594816*^9}, 3.636312159068121*^9, {3.636312331248557*^9, 
   3.6363124255238953`*^9}, 3.6363855129452477`*^9, {3.636391216454431*^9, 
   3.636391467489355*^9}, {3.6363918585692387`*^9, 3.63639186043509*^9}, {
   3.63639199101974*^9, 3.636392007043878*^9}, {3.636392427507909*^9, 
   3.636392430685143*^9}, {3.636392492020487*^9, 3.636392498284109*^9}, 
   3.6363925289309187`*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell["Quadrilateral parametric space", "Section",
 CellChangeTimes->{{3.636142521436479*^9, 3.636142524139866*^9}, {
  3.636143971451137*^9, 3.6361439829008636`*^9}}],

Cell[BoxData[{
 RowBox[{
  RowBox[{"If", "[", 
   RowBox[{"UseExactIntegration", ",", "\[IndentingNewLine]", 
    RowBox[{
     RowBox[{"TwoDimensionalQuadrilateral", "=", 
      RowBox[{"Table", "[", 
       RowBox[{
        RowBox[{"Integrate", "[", 
         RowBox[{
          RowBox[{"Total", "[", 
           RowBox[{"Flatten", "[", 
            RowBox[{"Table", "[", 
             RowBox[{
              RowBox[{"twodimensional", "/.", 
               RowBox[{"P", "\[Rule]", "p"}]}], ",", 
              RowBox[{"{", 
               RowBox[{"i", ",", "0", ",", "p"}], "}"}], ",", 
              RowBox[{"{", 
               RowBox[{"j", ",", "0", ",", "p"}], "}"}]}], "]"}], "]"}], 
           "]"}], ",", 
          RowBox[{"{", 
           RowBox[{"\[Xi]", ",", 
            RowBox[{"-", "1"}], ",", "1"}], "}"}], ",", 
          RowBox[{"{", 
           RowBox[{"\[Eta]", ",", 
            RowBox[{"-", "1"}], ",", "1"}], "}"}]}], "]"}], ",", 
        RowBox[{"{", 
         RowBox[{"p", ",", "0", ",", "Porder"}], "}"}]}], "]"}]}], ";"}], 
    "\[IndentingNewLine]", ",", "\[IndentingNewLine]", 
    RowBox[{
     RowBox[{"TwoDimensionalQuadrilateral", "=", 
      RowBox[{"Table", "[", 
       RowBox[{
        RowBox[{"NIntegrate", "[", 
         RowBox[{
          RowBox[{"Total", "[", 
           RowBox[{"Flatten", "[", 
            RowBox[{"Table", "[", 
             RowBox[{
              RowBox[{"twodimensional", "/.", 
               RowBox[{"P", "\[Rule]", "p"}]}], ",", 
              RowBox[{"{", 
               RowBox[{"i", ",", "0", ",", "p"}], "}"}], ",", 
              RowBox[{"{", 
               RowBox[{"j", ",", "0", ",", "p"}], "}"}]}], "]"}], "]"}], 
           "]"}], ",", 
          RowBox[{"{", 
           RowBox[{"\[Xi]", ",", 
            RowBox[{"-", "1"}], ",", "1"}], "}"}], ",", 
          RowBox[{"{", 
           RowBox[{"\[Eta]", ",", 
            RowBox[{"-", "1"}], ",", "1"}], "}"}], ",", 
          RowBox[{"Method", "\[Rule]", "method"}], ",", 
          RowBox[{"WorkingPrecision", "\[Rule]", "workpre"}]}], "]"}], ",", 
        RowBox[{"{", 
         RowBox[{"p", ",", "0", ",", "Porder"}], "}"}]}], "]"}]}], ";"}]}], 
   "\[IndentingNewLine]", "]"}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"file2DQ", "=", 
   RowBox[{"OpenWrite", "[", "\"\<Quadrilateral.txt\>\"", "]"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"Do", "[", " ", 
   RowBox[{
    RowBox[{"Write", "[", 
     RowBox[{"file2DQ", ",", 
      RowBox[{"TwoDimensionalQuadrilateral", "[", 
       RowBox[{"[", "n", "]"}], "]"}]}], "]"}], ",", 
    RowBox[{"{", 
     RowBox[{"n", ",", "1", ",", "Porder"}], "}"}]}], "]"}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"Close", "[", "file2DQ", "]"}], ";"}]}], "Input",
 CellChangeTimes->{{3.636143710281425*^9, 3.636143796641883*^9}, {
   3.6361438407822227`*^9, 3.636143866794297*^9}, {3.636143903500929*^9, 
   3.636143966659725*^9}, {3.63614407559863*^9, 3.6361440809457493`*^9}, 
   3.6361441283929443`*^9, {3.636144297998473*^9, 3.636144298326858*^9}, {
   3.636144377228163*^9, 3.636144447030261*^9}, {3.636144530736697*^9, 
   3.6361445355153418`*^9}, {3.636144576607875*^9, 3.6361445806709414`*^9}, 
   3.636145916325968*^9, {3.63614600328793*^9, 3.63614604833706*^9}, {
   3.63614610170501*^9, 3.63614611942173*^9}, 3.636146195006885*^9, {
   3.636146306762343*^9, 3.636146312729094*^9}, {3.636148211369488*^9, 
   3.636148223159886*^9}, 3.636155613618444*^9, 3.636294996902823*^9, {
   3.6362953223444567`*^9, 3.6362953679810038`*^9}, {3.6362955109890623`*^9, 
   3.6362955129866953`*^9}, {3.636295649937957*^9, 3.6362956531117907`*^9}, 
   3.6363121549181633`*^9, {3.636312482813697*^9, 3.636312490077046*^9}, {
   3.636393133473248*^9, 3.6363931347696867`*^9}}]
}, Open  ]],

Cell[CellGroupData[{

Cell["Tetrahedral parametric space", "Section",
 CellChangeTimes->{{3.636142521436479*^9, 3.636142524139866*^9}, {
  3.636143971451137*^9, 3.6361439829008636`*^9}, {3.636144893949394*^9, 
  3.636144894524077*^9}, {3.63614500985713*^9, 3.636145012168703*^9}}],

Cell[BoxData[{
 RowBox[{
  RowBox[{"If", "[", 
   RowBox[{"UseExactIntegration", ",", "\[IndentingNewLine]", 
    RowBox[{
     RowBox[{"ThreeDimensionalTetrahedral", "=", 
      RowBox[{"Table", "[", 
       RowBox[{
        RowBox[{"Integrate", "[", 
         RowBox[{
          RowBox[{"Total", "[", 
           RowBox[{"Flatten", "[", 
            RowBox[{"Table", "[", 
             RowBox[{
              RowBox[{"threedimensional", "/.", 
               RowBox[{"P", "\[Rule]", "p"}]}], ",", 
              RowBox[{"{", 
               RowBox[{"i", ",", "0", ",", "p"}], "}"}], ",", 
              RowBox[{"{", 
               RowBox[{"j", ",", "0", ",", "p"}], "}"}], ",", 
              RowBox[{"{", 
               RowBox[{"k", ",", "0", ",", "p"}], "}"}]}], "]"}], "]"}], 
           "]"}], ",", 
          RowBox[{"{", 
           RowBox[{"\[Xi]", ",", "0", ",", "1"}], "}"}], ",", 
          RowBox[{"{", 
           RowBox[{"\[Eta]", ",", "0", ",", 
            RowBox[{"1", "-", "\[Xi]"}]}], "}"}], ",", 
          RowBox[{"{", 
           RowBox[{"\[Zeta]", ",", "0", ",", 
            RowBox[{"1", "-", 
             RowBox[{"(", 
              RowBox[{"\[Xi]", "+", "\[Eta]"}], ")"}]}]}], "}"}]}], "]"}], 
        ",", 
        RowBox[{"{", 
         RowBox[{"p", ",", "0", ",", "Porder"}], "}"}]}], "]"}]}], ";"}], 
    "\[IndentingNewLine]", ",", "\[IndentingNewLine]", 
    RowBox[{
     RowBox[{"ThreeDimensionalTetrahedral", "=", 
      RowBox[{"Table", "[", 
       RowBox[{
        RowBox[{"NIntegrate", "[", 
         RowBox[{
          RowBox[{"Total", "[", 
           RowBox[{"Flatten", "[", 
            RowBox[{"Table", "[", 
             RowBox[{
              RowBox[{"threedimensional", "/.", 
               RowBox[{"P", "\[Rule]", "p"}]}], ",", 
              RowBox[{"{", 
               RowBox[{"i", ",", "0", ",", "p"}], "}"}], ",", 
              RowBox[{"{", 
               RowBox[{"j", ",", "0", ",", "p"}], "}"}], ",", 
              RowBox[{"{", 
               RowBox[{"k", ",", "0", ",", "p"}], "}"}]}], "]"}], "]"}], 
           "]"}], ",", 
          RowBox[{"{", 
           RowBox[{"\[Xi]", ",", "0", ",", "1"}], "}"}], ",", 
          RowBox[{"{", 
           RowBox[{"\[Eta]", ",", "0", ",", 
            RowBox[{"1", "-", "\[Xi]"}]}], "}"}], ",", 
          RowBox[{"{", 
           RowBox[{"\[Zeta]", ",", "0", ",", 
            RowBox[{"1", "-", 
             RowBox[{"(", 
              RowBox[{"\[Xi]", "+", "\[Eta]"}], ")"}]}]}], "}"}], ",", 
          RowBox[{"Method", "\[Rule]", "method"}], ",", 
          RowBox[{"WorkingPrecision", "\[Rule]", "workpre"}]}], "]"}], ",", 
        RowBox[{"{", 
         RowBox[{"p", ",", "0", ",", "Porder"}], "}"}]}], "]"}]}], ";"}]}], 
   "\[IndentingNewLine]", "]"}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"file3DT", "=", 
   RowBox[{"OpenWrite", "[", "\"\<Tetrahedral.txt\>\"", "]"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"Do", "[", " ", 
   RowBox[{
    RowBox[{"Write", "[", 
     RowBox[{"file3DT", ",", 
      RowBox[{"ThreeDimensionalTetrahedral", "[", 
       RowBox[{"[", "n", "]"}], "]"}]}], "]"}], ",", 
    RowBox[{"{", 
     RowBox[{"n", ",", "1", ",", "Porder"}], "}"}]}], "]"}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"Close", "[", "file3DT", "]"}], ";"}]}], "Input",
 CellChangeTimes->{{3.636143710281425*^9, 3.636143796641883*^9}, {
   3.6361438407822227`*^9, 3.636143866794297*^9}, {3.636143903500929*^9, 
   3.636143966659725*^9}, {3.63614407559863*^9, 3.6361440809457493`*^9}, 
   3.6361441283929443`*^9, {3.636144297998473*^9, 3.636144298326858*^9}, {
   3.636144377228163*^9, 3.636144447030261*^9}, {3.636144530736697*^9, 
   3.6361445355153418`*^9}, {3.636144576607875*^9, 3.6361445806709414`*^9}, {
   3.636144913816265*^9, 3.636144983543839*^9}, {3.636145016146476*^9, 
   3.636145095138297*^9}, {3.636145920924781*^9, 3.63614592217844*^9}, {
   3.63614607944349*^9, 3.636146080156755*^9}, {3.636146123066762*^9, 
   3.636146125956128*^9}, 3.636146196458868*^9, 3.636146317923321*^9, {
   3.636148211374632*^9, 3.636148223164316*^9}, {3.636155246561912*^9, 
   3.63615524880228*^9}, {3.636155416614307*^9, 3.636155452702405*^9}, {
   3.636155485464765*^9, 3.636155485774501*^9}, {3.636155582714467*^9, 
   3.636155583939863*^9}, 3.6361556176891108`*^9, {3.636295440624683*^9, 
   3.63629548676888*^9}, {3.6362956451898212`*^9, 3.636295703637721*^9}, 
   3.636312598049383*^9, {3.636393136511153*^9, 3.636393137554449*^9}}]
}, Open  ]],

Cell[CellGroupData[{

Cell["Cube parametric space", "Section",
 CellChangeTimes->{{3.636142521436479*^9, 3.636142524139866*^9}, {
  3.636143971451137*^9, 3.6361439829008636`*^9}, {3.636144893949394*^9, 
  3.636144894524077*^9}}],

Cell[BoxData[{
 RowBox[{
  RowBox[{"If", "[", 
   RowBox[{"UseExactIntegration", ",", "\[IndentingNewLine]", 
    RowBox[{
     RowBox[{"ThreeDimensionalCube", "=", 
      RowBox[{"Table", "[", 
       RowBox[{
        RowBox[{"Integrate", "[", 
         RowBox[{
          RowBox[{"Total", "[", 
           RowBox[{"Flatten", "[", 
            RowBox[{"Table", "[", 
             RowBox[{
              RowBox[{"threedimensional", "/.", 
               RowBox[{"P", "\[Rule]", "p"}]}], ",", 
              RowBox[{"{", 
               RowBox[{"i", ",", "0", ",", "p"}], "}"}], ",", 
              RowBox[{"{", 
               RowBox[{"j", ",", "0", ",", "p"}], "}"}], ",", 
              RowBox[{"{", 
               RowBox[{"k", ",", "0", ",", "p"}], "}"}]}], "]"}], "]"}], 
           "]"}], ",", 
          RowBox[{"{", 
           RowBox[{"\[Xi]", ",", 
            RowBox[{"-", "1"}], ",", "1"}], "}"}], ",", 
          RowBox[{"{", 
           RowBox[{"\[Eta]", ",", 
            RowBox[{"-", "1"}], ",", "1"}], "}"}], ",", 
          RowBox[{"{", 
           RowBox[{"\[Zeta]", ",", 
            RowBox[{"-", "1"}], ",", "1"}], "}"}]}], "]"}], ",", 
        RowBox[{"{", 
         RowBox[{"p", ",", "0", ",", "Porder"}], "}"}]}], "]"}]}], ";"}], 
    "\[IndentingNewLine]", ",", "\[IndentingNewLine]", 
    RowBox[{
     RowBox[{"ThreeDimensionalCube", "=", 
      RowBox[{"Table", "[", 
       RowBox[{
        RowBox[{"NIntegrate", "[", 
         RowBox[{
          RowBox[{"Total", "[", 
           RowBox[{"Flatten", "[", 
            RowBox[{"Table", "[", 
             RowBox[{
              RowBox[{"threedimensional", "/.", 
               RowBox[{"P", "\[Rule]", "p"}]}], ",", 
              RowBox[{"{", 
               RowBox[{"i", ",", "0", ",", "p"}], "}"}], ",", 
              RowBox[{"{", 
               RowBox[{"j", ",", "0", ",", "p"}], "}"}], ",", 
              RowBox[{"{", 
               RowBox[{"k", ",", "0", ",", "p"}], "}"}]}], "]"}], "]"}], 
           "]"}], ",", 
          RowBox[{"{", 
           RowBox[{"\[Xi]", ",", 
            RowBox[{"-", "1"}], ",", "1"}], "}"}], ",", 
          RowBox[{"{", 
           RowBox[{"\[Eta]", ",", 
            RowBox[{"-", "1"}], ",", "1"}], "}"}], ",", 
          RowBox[{"{", 
           RowBox[{"\[Zeta]", ",", 
            RowBox[{"-", "1"}], ",", "1"}], "}"}], ",", 
          RowBox[{"Method", "\[Rule]", "method"}], ",", 
          RowBox[{"WorkingPrecision", "\[Rule]", "workpre"}]}], "]"}], ",", 
        RowBox[{"{", 
         RowBox[{"p", ",", "0", ",", "Porder"}], "}"}]}], "]"}]}], ";"}]}], 
   "\[IndentingNewLine]", "]"}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"file3DC", "=", 
   RowBox[{"OpenWrite", "[", "\"\<Cube.txt\>\"", "]"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"Do", "[", " ", 
   RowBox[{
    RowBox[{"Write", "[", 
     RowBox[{"file3DC", ",", 
      RowBox[{"ThreeDimensionalCube", "[", 
       RowBox[{"[", "n", "]"}], "]"}]}], "]"}], ",", 
    RowBox[{"{", 
     RowBox[{"n", ",", "1", ",", "Porder"}], "}"}]}], "]"}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"Close", "[", "file3DC", "]"}], ";"}]}], "Input",
 CellChangeTimes->{{3.636143710281425*^9, 3.636143796641883*^9}, {
   3.6361438407822227`*^9, 3.636143866794297*^9}, {3.636143903500929*^9, 
   3.636143966659725*^9}, {3.63614407559863*^9, 3.6361440809457493`*^9}, 
   3.6361441283929443`*^9, {3.636144297998473*^9, 3.636144298326858*^9}, {
   3.636144377228163*^9, 3.636144447030261*^9}, {3.636144530736697*^9, 
   3.6361445355153418`*^9}, {3.636144576607875*^9, 3.6361445806709414`*^9}, {
   3.636144913816265*^9, 3.636144983543839*^9}, 3.636145926146817*^9, {
   3.636146127096863*^9, 3.63614612911866*^9}, 3.6361461978489122`*^9, 
   3.636146321530396*^9, {3.636148211380176*^9, 3.63614822316908*^9}, 
   3.6361556220132113`*^9, {3.6362957698727827`*^9, 3.63629580202559*^9}, 
   3.6363126360193653`*^9, {3.636393139504177*^9, 3.636393141202828*^9}}]
}, Open  ]],

Cell[CellGroupData[{

Cell["Pyramid parametric space", "Section",
 CellChangeTimes->{{3.636142521436479*^9, 3.636142524139866*^9}, {
  3.636143971451137*^9, 3.6361439829008636`*^9}, {3.636144893949394*^9, 
  3.636144894524077*^9}, {3.636145161691929*^9, 3.6361451634798594`*^9}}],

Cell[BoxData[{
 RowBox[{
  RowBox[{"If", "[", 
   RowBox[{"UseExactIntegration", ",", "\[IndentingNewLine]", 
    RowBox[{
     RowBox[{"ThreeDimensionalPyramid", "=", 
      RowBox[{"Table", "[", 
       RowBox[{
        RowBox[{"Integrate", "[", 
         RowBox[{
          RowBox[{"Total", "[", 
           RowBox[{"Flatten", "[", 
            RowBox[{"Table", "[", 
             RowBox[{
              RowBox[{"threedimensional", "/.", 
               RowBox[{"P", "\[Rule]", "p"}]}], ",", 
              RowBox[{"{", 
               RowBox[{"i", ",", "0", ",", "p"}], "}"}], ",", 
              RowBox[{"{", 
               RowBox[{"j", ",", "0", ",", "p"}], "}"}], ",", 
              RowBox[{"{", 
               RowBox[{"k", ",", "0", ",", "p"}], "}"}]}], "]"}], "]"}], 
           "]"}], ",", 
          RowBox[{"{", 
           RowBox[{"\[Zeta]", ",", "0", ",", "1"}], "}"}], ",", 
          RowBox[{"{", 
           RowBox[{"\[Xi]", ",", 
            RowBox[{"\[Zeta]", "-", "1"}], ",", 
            RowBox[{"1", "-", "\[Zeta]"}]}], "}"}], ",", 
          RowBox[{"{", 
           RowBox[{"\[Eta]", ",", 
            RowBox[{"\[Zeta]", "-", "1"}], ",", 
            RowBox[{"1", "-", "\[Zeta]"}]}], "}"}]}], "]"}], ",", 
        RowBox[{"{", 
         RowBox[{"p", ",", "0", ",", "Porder"}], "}"}]}], "]"}]}], ";"}], 
    "\[IndentingNewLine]", ",", "\[IndentingNewLine]", 
    RowBox[{
     RowBox[{"ThreeDimensionalPyramid", "=", 
      RowBox[{"Table", "[", 
       RowBox[{
        RowBox[{"NIntegrate", "[", 
         RowBox[{
          RowBox[{"Total", "[", 
           RowBox[{"Flatten", "[", 
            RowBox[{"Table", "[", 
             RowBox[{
              RowBox[{"threedimensional", "/.", 
               RowBox[{"P", "\[Rule]", "p"}]}], ",", 
              RowBox[{"{", 
               RowBox[{"i", ",", "0", ",", "p"}], "}"}], ",", 
              RowBox[{"{", 
               RowBox[{"j", ",", "0", ",", "p"}], "}"}], ",", 
              RowBox[{"{", 
               RowBox[{"k", ",", "0", ",", "p"}], "}"}]}], "]"}], "]"}], 
           "]"}], ",", 
          RowBox[{"{", 
           RowBox[{"\[Zeta]", ",", "0", ",", "1"}], "}"}], ",", 
          RowBox[{"{", 
           RowBox[{"\[Xi]", ",", 
            RowBox[{"\[Zeta]", "-", "1"}], ",", 
            RowBox[{"1", "-", "\[Zeta]"}]}], "}"}], ",", 
          RowBox[{"{", 
           RowBox[{"\[Eta]", ",", 
            RowBox[{"\[Zeta]", "-", "1"}], ",", 
            RowBox[{"1", "-", "\[Zeta]"}]}], "}"}], ",", 
          RowBox[{"Method", "\[Rule]", "method"}], ",", 
          RowBox[{"WorkingPrecision", "\[Rule]", "workpre"}]}], "]"}], ",", 
        RowBox[{"{", 
         RowBox[{"p", ",", "0", ",", "Porder"}], "}"}]}], "]"}]}], ";"}]}], 
   "\[IndentingNewLine]", "]"}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"file3DP", "=", 
   RowBox[{"OpenWrite", "[", "\"\<Pyramid.txt\>\"", "]"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"Do", "[", " ", 
   RowBox[{
    RowBox[{"Write", "[", 
     RowBox[{"file3DP", ",", 
      RowBox[{"ThreeDimensionalPyramid", "[", 
       RowBox[{"[", "n", "]"}], "]"}]}], "]"}], ",", 
    RowBox[{"{", 
     RowBox[{"n", ",", "1", ",", "Porder"}], "}"}]}], "]"}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"Close", "[", "file3DP", "]"}], ";"}]}], "Input",
 CellChangeTimes->{{3.636143710281425*^9, 3.636143796641883*^9}, {
   3.6361438407822227`*^9, 3.636143866794297*^9}, {3.636143903500929*^9, 
   3.636143966659725*^9}, {3.63614407559863*^9, 3.6361440809457493`*^9}, 
   3.6361441283929443`*^9, {3.636144297998473*^9, 3.636144298326858*^9}, {
   3.636144377228163*^9, 3.636144447030261*^9}, {3.636144530736697*^9, 
   3.6361445355153418`*^9}, {3.636144576607875*^9, 3.6361445806709414`*^9}, {
   3.636144913816265*^9, 3.636144983543839*^9}, {3.63614516799726*^9, 
   3.63614523143382*^9}, 3.636145265661599*^9, 3.636145928983652*^9, {
   3.636146131136647*^9, 3.636146133022647*^9}, 3.636146199681245*^9, 
   3.636146324042652*^9, {3.6361482113864*^9, 3.636148223174758*^9}, 
   3.636155624787024*^9, {3.6362958303439713`*^9, 3.636295851984017*^9}, 
   3.6363126394641857`*^9, {3.636393145760337*^9, 3.636393147091755*^9}}]
}, Open  ]],

Cell[CellGroupData[{

Cell["Prism parametric space", "Section",
 CellChangeTimes->{{3.636142521436479*^9, 3.636142524139866*^9}, {
  3.636143971451137*^9, 3.6361439829008636`*^9}, {3.636144893949394*^9, 
  3.636144894524077*^9}, {3.636145161691929*^9, 3.6361451634798594`*^9}, {
  3.636145276828629*^9, 3.6361452805809383`*^9}}],

Cell[BoxData[{
 RowBox[{
  RowBox[{"If", "[", 
   RowBox[{"UseExactIntegration", ",", "\[IndentingNewLine]", 
    RowBox[{
     RowBox[{"ThreeDimensionalPrism", "=", 
      RowBox[{"Table", "[", 
       RowBox[{
        RowBox[{"Integrate", "[", 
         RowBox[{
          RowBox[{"Total", "[", 
           RowBox[{"Flatten", "[", 
            RowBox[{"Table", "[", 
             RowBox[{
              RowBox[{"threedimensional", "/.", 
               RowBox[{"P", "\[Rule]", "p"}]}], ",", 
              RowBox[{"{", 
               RowBox[{"i", ",", "0", ",", "p"}], "}"}], ",", 
              RowBox[{"{", 
               RowBox[{"j", ",", "0", ",", "p"}], "}"}], ",", 
              RowBox[{"{", 
               RowBox[{"k", ",", "0", ",", "p"}], "}"}]}], "]"}], "]"}], 
           "]"}], ",", 
          RowBox[{"{", 
           RowBox[{"\[Xi]", ",", "0", ",", "1"}], "}"}], ",", 
          RowBox[{"{", 
           RowBox[{"\[Eta]", ",", "0", ",", 
            RowBox[{"1", "-", "\[Xi]"}]}], "}"}], ",", 
          RowBox[{"{", 
           RowBox[{"\[Zeta]", ",", 
            RowBox[{"-", "1"}], ",", "1"}], "}"}]}], "]"}], ",", 
        RowBox[{"{", 
         RowBox[{"p", ",", "0", ",", "Porder"}], "}"}]}], "]"}]}], ";"}], 
    "\[IndentingNewLine]", ",", "\[IndentingNewLine]", 
    RowBox[{
     RowBox[{"ThreeDimensionalPrism", "=", 
      RowBox[{"Table", "[", 
       RowBox[{
        RowBox[{"NIntegrate", "[", 
         RowBox[{
          RowBox[{"Total", "[", 
           RowBox[{"Flatten", "[", 
            RowBox[{"Table", "[", 
             RowBox[{
              RowBox[{"threedimensional", "/.", 
               RowBox[{"P", "\[Rule]", "p"}]}], ",", 
              RowBox[{"{", 
               RowBox[{"i", ",", "0", ",", "p"}], "}"}], ",", 
              RowBox[{"{", 
               RowBox[{"j", ",", "0", ",", "p"}], "}"}], ",", 
              RowBox[{"{", 
               RowBox[{"k", ",", "0", ",", "p"}], "}"}]}], "]"}], "]"}], 
           "]"}], ",", 
          RowBox[{"{", 
           RowBox[{"\[Xi]", ",", "0", ",", "1"}], "}"}], ",", 
          RowBox[{"{", 
           RowBox[{"\[Eta]", ",", "0", ",", 
            RowBox[{"1", "-", "\[Xi]"}]}], "}"}], ",", 
          RowBox[{"{", 
           RowBox[{"\[Zeta]", ",", 
            RowBox[{"-", "1"}], ",", "1"}], "}"}], ",", 
          RowBox[{"Method", "\[Rule]", "method"}], ",", 
          RowBox[{"WorkingPrecision", "\[Rule]", "workpre"}]}], "]"}], ",", 
        RowBox[{"{", 
         RowBox[{"p", ",", "0", ",", "Porder"}], "}"}]}], "]"}]}], ";"}]}], 
   "\[IndentingNewLine]", "]"}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"file3DPr", "=", 
   RowBox[{"OpenWrite", "[", "\"\<Prism.txt\>\"", "]"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"Do", "[", " ", 
   RowBox[{
    RowBox[{"Write", "[", 
     RowBox[{"file3DPr", ",", 
      RowBox[{"ThreeDimensionalPrism", "[", 
       RowBox[{"[", "n", "]"}], "]"}]}], "]"}], ",", 
    RowBox[{"{", 
     RowBox[{"n", ",", "1", ",", "Porder"}], "}"}]}], "]"}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"Close", "[", "file3DPr", "]"}], ";"}]}], "Input",
 CellChangeTimes->{{3.636143710281425*^9, 3.636143796641883*^9}, {
   3.6361438407822227`*^9, 3.636143866794297*^9}, {3.636143903500929*^9, 
   3.636143966659725*^9}, {3.63614407559863*^9, 3.6361440809457493`*^9}, 
   3.6361441283929443`*^9, {3.636144297998473*^9, 3.636144298326858*^9}, {
   3.636144377228163*^9, 3.636144447030261*^9}, {3.636144530736697*^9, 
   3.6361445355153418`*^9}, {3.636144576607875*^9, 3.6361445806709414`*^9}, {
   3.636144913816265*^9, 3.636144983543839*^9}, {3.63614516799726*^9, 
   3.63614523143382*^9}, {3.636145265661599*^9, 3.6361453349005547`*^9}, 
   3.6361456565107*^9, 3.636145933512272*^9, {3.636146135385109*^9, 
   3.636146137239118*^9}, 3.636146201008932*^9, 3.6361463279513817`*^9, {
   3.636148211393833*^9, 3.636148223179514*^9}, 3.636155627442904*^9, {
   3.6362958699815893`*^9, 3.636295899757338*^9}, 3.636312643366768*^9, {
   3.636393148470276*^9, 3.636393149723948*^9}}]
}, Open  ]]
}, Open  ]]
},
WindowSize->{2102, 869},
WindowMargins->{{Automatic, 45}, {Automatic, 89}},
FrontEndVersion->"10.0 for Mac OS X x86 (32-bit, 64-bit Kernel) (September 9, \
2014)",
StyleDefinitions->"Default.nb"
]
(* End of Notebook Content *)

(* Internal cache information *)
(*CellTagsOutline
CellTagsIndex->{}
*)
(*CellTagsIndex
CellTagsIndex->{}
*)
(*NotebookFileOutline
Notebook[{
Cell[CellGroupData[{
Cell[580, 22, 107, 1, 92, "Title"],
Cell[690, 25, 3002, 50, 97, "Input"],
Cell[CellGroupData[{
Cell[3717, 79, 96, 1, 64, "Section"],
Cell[3816, 82, 2136, 62, 134, "Input"]
}, Open  ]],
Cell[CellGroupData[{
Cell[5989, 149, 205, 3, 64, "Section"],
Cell[6197, 154, 4238, 88, 131, "Input"]
}, Open  ]],
Cell[CellGroupData[{
Cell[10472, 247, 258, 3, 64, "Section"],
Cell[10733, 252, 4278, 93, 148, "Input"]
}, Open  ]],
Cell[CellGroupData[{
Cell[15048, 350, 166, 2, 64, "Section"],
Cell[15217, 354, 3778, 87, 148, "Input"]
}, Open  ]],
Cell[CellGroupData[{
Cell[19032, 446, 258, 3, 64, "Section"],
Cell[19293, 451, 4493, 102, 148, "Input"]
}, Open  ]],
Cell[CellGroupData[{
Cell[23823, 558, 206, 3, 64, "Section"],
Cell[24032, 563, 3967, 94, 148, "Input"]
}, Open  ]],
Cell[CellGroupData[{
Cell[28036, 662, 257, 3, 64, "Section"],
Cell[28296, 667, 4190, 97, 148, "Input"]
}, Open  ]],
Cell[CellGroupData[{
Cell[32523, 769, 306, 4, 64, "Section"],
Cell[32832, 775, 4055, 94, 206, "Input"]
}, Open  ]]
}, Open  ]]
}
]
*)

(* End of internal cache information *)
